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1. Introduction 



After the breakthroughs of 2005 [U EJ [3] with the fully non-linear dynamical numerical 
simulation of the inspiral, merger and ringdown of black- hole binaries (BHBs), there 
were many important advances in the understanding of black-hole physics. Indeed, the 
discovery and modeling of very large recoil velocities [H [5] acquired by the final remnant 
of spinning BHB mergers have attracted a lot of interest among astrophysicists. 

Empirical formulae for the final remnant black hole recoil velocity (also mass and 
spin) from merging black- hole binaries were obtained in [6] (and references therein), 
where post-Newtonian results (TJ Ej were used as guide to model the recoil dependence 
on the physical parameters of the progenitor BHB [6]. On the other hand, there is also 
a long history of recoil studies in black hole perturbation theory (e.g., see [10]), and 
recently, an analytic treatment of the linear momentum flux of the plunge of a particle 
into a Kerr black hole has been considered in 

In [12], we extended the Regge-Wheeler-Zerilli (RWZ) equations [131 Ej f° r the 
Schwarzschild perturbations by including, perturbatively, a term linear in the spin of 
the larger black hole (SRWZ formulation). We have found good agreement in the full 
numerical and perturbative waveforms for intermediate mass ratio black-hole binaries, 
reaching 99.5% matching for the leading (£,m) = (2,2) mode [151 E2] and recently 
simulated a mass ratio 100 : 1 BHB merger [16]. This formalism can be considered as 
a extension of [17] in the close limit. Here, using the slow motion approximation, we 
derive the evolution of the linear momentum for binary systems in the SRWZ formalism 
analytically and compare to the post-Newtonian (PN) expansions. 



2. Formulation 



2.1. Spin as a perturbation 

Our goal is to analytically model the waveforms from a particle with mass fi orbiting 
around a spinning black hole with mass M (where [i -C M). In [12], we considered the 
Kerr metric up to 0(a l ), where a denotes the spin of the BH which has the dimension 
of mass, and the spin direction is along the z-axis. In this paper, we set the spin along 
the x-axis, and the metric is given by 

ds 2 - 1 1 \M2 , li \ 



\ d t 2 + h - filij dr 2 + r 2 (d9 2 + sin 2 

dt (sin (J)d9 + sin 9 cos 9 cos (j)d(j)) + 0(a 2 ) . (1) 



r 

in the Boyer-Lindquist radial coordinate. The last term in the right hand side of the 
above equation is treated as a perturbation in the Schwarzschild black hole background. 

n - n Sch 4- h^' spin ^ ■ 

yfiv y^u ~ "'tin i 

no o C 

,(l,spin) _ r(l,spm) _ f>fx . , ,(l,spin) _ ,(l,spin) _ „• a „„„ a nr .„ A. /o\ 

h t9 - h et sin , hfy - tifa sin 9 cos 9cos(p, (2) 
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where S x = Ma. By using the expansion defined by the tensor harmonics given in 
we find that the coefficients of the tensor harmonics are given by 




(l,spin) 
01-1 



(^ r ) = \hr — ■ 



(3) 



2.2. SEW Z formulation 

We treat the coupling between the spin discussed above and first order metric 
perturbation (gravitational radiation) as a second order perturbation calculation. The 
Einstein equations up to the second perturbative order are formally written as 

&$\&\ + G^ (4) 

where the energy- momentum tensor T^ u = T$ + T^}, and hfy and hffi are the first 
and second order metric perturbations, respectively. We consider that the second 
order metric perturbation, /i( 2 > wave ) is created by the spin h^' sipm ' and radiation /?,( 1 ' wavc ) 
couplings. In this case, we solve 

G (l)^(l,wave) ]=87rT ^ ? (5) 
G (l) ^(2,wave)] = _ q{2) ^(l.wave) ? ^(l,spin)] ? (q) 

up to 0(a l ). Here the square of the first order radiation has been ignored. We solve the 
above equations in the tensor harmonics expansion as in the RWZ formalism. 



2.3. Gravitational radiation and linear momentum evolution 



Gravitational wave modes in the RWZ formalism are given by the metric perturbation 
under an asymptotic flat (AF) gauge 



£m I 



-£(£+!)(£ -!)(£ + 2) 



1/2 



G 



(n)AF 



Cm 



(in 



[2£(£+l)(£-l)(£+2)] 1/2 (n)AF 



2r s 



2£m 



\t,r)d 



I'm 



(7) 



where superscripts (n) (n = 1, 2) denote the perturbative order, and f£ m and di m 
are tensor harmonics. The even parity mode G^ AF and odd parity mode li^j^f are 

MCeven) _ A n/ ,(n)(odd) as 



expressed by the Regge-Wheeler-Zerilli functions ip y ^ yavaLL ' and ijj em 



G 



(n)AF 
I'm 



(t,r) 



(n)(even) 



(t,r) 



h 



(n)AF 
2£m 



(t,r) 



irip^ {odd \t,r) 



(8) 



(See the definition of the tensor harmonics expansion given in |18j . and [12] for more 
precise discussion of the second perturbative order.) 

Next, we discuss the time evolutions of the linear momentum of the binary system. 
This is obtained from the following expression (e.g. see Eq. (2.17) in jlO]). 

1 



Pi= - 



327T 



dQr 2 ni (h afS . t h 



TT 



(9) 
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where rij = (sin 9 cos 0, sin 9 sin 0, COS6 1 ). The subscript TT means the transverse 
traceless part and we may use the metric in Eq. ([7]) under the AF gauge. Then, the 
above equation is written as 

£m I'm' 



r 2 Gt )AF (t, r) G^ F (f r) - 1 (*, /^ft r) ) Pf («W) 

(Gt )AF (t, (* , r) - h^ F (t, r) \f r)) P, C («W) 



(10) 



where 



Pf{imi'm') = fain ni (w em W e , m , + — L-X £m Jfy, 
J V sm 

P. c (Wm') = fdnn,-^— (W em X e , m , - X tm W em >) . (11) 
j sm v 

We note that Pf vanishes for £ — £' = 2k (k: integer) because of the parity of the 
integration, while Pf vanishes for £ — £' = Ik + 1 (k: integer). 

In the following, we discuss only the leading order contribution of gravitational 
waves to the linear momentum evolution in the slow motion approximation. In practice, 
the combination of the first order £ = 2, m = and £ = 3, m = even parity modes 
produces the leading order contribution for the non-spinning case, and the combination 
of £ = 2, m = even parity and the £ = 2, m = ±1 odd parity modes becomes the 
dominant contribution of spin in the situation discussed below. 



3. Evolution of the linear momentum 



3.1. Point particle's motion 



We consider a particle falling radially into a Schwarzschild black hole as the first 
order source. Assuming O(r) = 3>(t) = [p, where the particle's location is given 
by {T(t), R(t), 6(r), $(r)}, the equation of motion of the particle is 

(olR\ 2 ( 2M\ 3 1 / 2M\ 2 



V*J = -\}—r) & + V-ir) • (12) 

where R = R(t) is the location of the particle and the energy E is written by 

/ 2M\ alT(r) , , 



We also use 



ol 2 R 3 / 2M\ 2 M ( 2M\ M 



ttJ w + 2 V—r)b>' (14) 

to simplify the calculations below. In the slow motion approximation where we consider 
alR/alt < 1 and M/R < 1, we have E — 1 and d 2 R/dt 2 = -M/R 2 . 

| Although we write $(r) = here, $(r) is arbitrary at 0(t) = 0. 



Perturbative effects of spinning black holes 



5 



The tensor harmonics coefficients of the first order stress-energy tensor become 



A$L(t, r) = n R ^ t) m 2M (r r 2 4 M)2 *(r - i2(t)) >7 m (0, 0) , 

41 (t, r) = i /i ^ <*(r - R(t)) F/ m (0, 0) . (15) 

Otherwise the coefficients are zero. Here, only the m = modes have a non-zero value. 

We note that the black-hole spin effect in the equations of motion of the particle 
does not contribute to leading order in the slow motion approximation. Therefore, we 
do not consider any effect of the background spin on the particle's trajectory in this 
paper. 

3.2. First order I = 2, m = and i = 3, m = 0, even parity perturbation 



In order to calculate the radiative even parity modes, we use the Zerilli equation. For 
example, the Zerilli equation for the I = 2, m = mode is given by 

4o )(even) (*,r) 



d 2 d 2 (r-2M)(4r 3 + 4r 2 M + 6rM 2 + 3M 3 ) 
6 



dt 2 dr* 2 r 4 (2r + 3M) 2 

4 7r u . „, ,, Air a d . . 

-8 (r - R(t)) + -—6 (r - R(t)) 



m - v 3 , r -v W). (16) 

In the right hand side of the above equation, we used the slow motion approximation, i.e., 
the characteristic orbital velocity v <C 1 where the particle's velocity dR(t)/dt ~ v and 
the potential M/R(t) ~ v 2 , and only consider leading-order terms in this approximation. 

The Zerilli function is obtained by the Green's function method. In practice, we 
used the Fourier transformation and prepared the Green's function in the frequency 
domain. We consider only the outside Green's function of the particle's location in the 
M — > limit. This is given by 

G(r, r') = iuj 2 {ur') h£\ur) 6(r - r') . (17) 

Here, we have asymptotic behaviors, ^(wr) = (wr) 2 /15 for small r, and h^\ur) = 
i exp(iu)r) / (cur) for large r. Finally, we obtain 

&™\t,r) = l 4^Ji? -^W o (0,0) , (18) 



15 r V R 
where the argument of R is the retarded time (t — r). 

In the same way, we calculate the t = 3, m = even parity mode. In the leading 
order of the slow motion approximation, the source term of the Zerilli equation becomes 



°30 



(l)(even) 



8 7T U r . n , ,, 47T LL d . - n , 

-5 (r - R(t)) + -—8 (r - R(t)) 



Y*(0,0). (19) 

Then, using the Green's function method with asymptotic behaviors, jz{ujr) = 
(wr) 3 /105 for small r, and h^\ur) = exp(iur) / (ur) for large r. we obtain 

^(even) ^ r) = 16W# _ ^ ^ (0> Q) . (2Q) 



105 ^ V R 
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3. 3. First order dipole, even parity perturbation 

Although the first order dipole {£ = 1) even parity mode in the vacuum regions can be 
completely eliminated in the center of mass coordinate system, the metric perturbation 
is not globally pure gauge [19] . The coupling between this mode and the black-hole spin 
creates the leading order spin effect on the linear momentum evolution. We therefore 
have to include the I = 1 mode contributions. 

For this mode, the metric perturbations are given by 

hff = (i_H^)^¥o(*,r)aoio-V^i^¥o(*^)aiio+(l-^) ^fiftO aio 

~ h { ^ e) (t, r) h Q10 + - r h$$\t, r) b 10 + V2Kg\t, r) g 10 . (21) 
The generators of the gauge transformation are by, 

x» -> x" + eii } f (x a ) ; 

& = {v o ( %,r)F 10 (M), ^ (1) (t,r)y 10 (M, 

vM(t,r)d e Y 10 (9,t), y«(t,r)^4^l (22) 

sin 6 J 

where Vq 1 ^, Vi and are three degrees of gauge freedom in the 1=1 mode. The 
metric perturbations transform under the above gauge transformation from a gauge (G) 
to a gauge (C) as 

H®? (*, = Sfii? (t, r) + 2 1 C°" G ' («, r) + 2 -^JL—V^ («, r) , 
ff'Kffcr) = ffffftr) + LZM l-v^ ( t , r ) - _L_ (t,r) , 

<» G '(t,r) = K« G ( f ,r) - V>°- G ' (t,r) + 2y« G - G ' (t,r) , 

r 

« |G '(*,r) = «> G (i,r) + ZZ™vP*M ( t , r ) - r ^<» G -°' (t,r) , 

r or 

^(*.r) = & )G (t,r) - _L_V' 1 <1)G "* G ' (t,r) -r^' 1 '^' (t,r) . (23) 

When we choose the gauge so that hoio Z = ^110 = ^io^ Z = 0; where the suffix 
Z stands for the Zerilli gauge [H], we obtain the metric perturbations 

H ™ {t > r) = 3M(r- E 2My + MWt) ~ 2M) ) 9{r ~ m) Y ^ 0) ' 

H ^)--J^^0(r-R( t )) ^(0,0), 

^(i, r) = ( /_^ )2 W) - 2M) 0(r - fl(t)) ^(0, 0) . (24) 
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We note that the metric perturbations are not under the AF gauge. For the above 
metric perturbation, if we consider the gauge transformation, 

V m^( t r ) - r3 dR ^ P(r - R(t)) Y* (0 0) 

V ° {t ' r} ~ 3M (r-2M) J dt ^ ^ y io( u ^' 

vl 1)Z ^ D (f r) = _ r 2M) (i2(t) - 2 M) 0(r - fl(t)) 3? (0, 0) , 

W 1)Z ^ D (t, r) = (r _ 1 2M) (i2(t) - 2 M) 0(r - i2(t)) 1^(0, 0) , (25) 

we obtain the singular metric perturbation at the particle's location. 

H W (tr) - 8 ^ E R W ( dR ®\ 2 5( r R(t))Y* (0 0) 
^oio V,r) - 3M {m _ 2M)2 [ dt ) *ir R{t))Y 10 {0,0), 

-W^^^^-^)^(o,o), 

H$?(t,r) = ^R(t)5(r-R(t))Y* (0,0), 
h^ )D (f r) = - ^R(t)^ 6(r - *(,)) F; o (0, 0) , 
h^ )D (t, r) = ^R(t? 6(r - R(t)) Y* Q (0, 0) , 

KT{t,r) =0, (26) 

The above means that although the metric perturbations vanish in the vacuum regions, 
there are some contributions at the location of the particle. This gauge choice has been 
discussed in [Tj5] for circular orbit. The coordinate system under this gauge condition 
can be considered as the center of mass system which is suitable for the analysis of 
the evolution of the linear momentum. Therefore, we use this gauge and its metric 
perturbations here. It is noted that we also have another gauge choice where the metric 
perturbations become C° at the particle's location [2D] . 

3.4- Second order I = 2, m = ±1, odd parity perturbation 

In the previous subsections, we focused on the first perturbative order. Here we treat 
the second perturbative order using the SRWZ formalism. We focus on the coupling 
between the first order £ — 1, m = even parity mode (with parity (— l) 1 ) and the spin 
effect of the central black hole in Eq. which is given by the £ — 1, m = ±1 odd parity 
modes (with parity (— 1) 1+1 ). This coupling creates the second order £ — 2, m — ±1 
odd parity perturbation (with parity (— 1) 2+1 ). 

The Regge- Wheeler function for this odd parity perturbation satisfies 



d 2 d 2 (r-2M)(r-M) 
di 2 ~ + dr l2 ~ 7 4 



^r\t,r)=S$r(t,r). (27) 



The source term is derived from the effective stress-energy tensor, 

T^ cff ) = - -L^g) [h^ ole \ h {1 ' sp[n) ] , (28) 
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by using the same tensor harmonics expansion as for the first perturbative order. When 
we consider the leading order in the slow motion approximation, the source term S^±[ 
becomes 

* (rW S <r ' m) ~ W)T/ (r * m) + J? (r - R{t)) ) • (29) 
if we use the first order £ = 1, m = even parity mode under the D gauge. From this 
source term, we obtain the Regge- Wheeler function 

*gr>(t,r)= ±i^§y.(0,0), (30) 

by using the same Green's function method discussed in Subsection 13.21 Here, we have 
written S 2 = S x . 

3.5. Spinning particle orbiting around a black hole 

In the above subsection, we considered the black hole with mass M that has a spin along 
the x-direction. Here, we introduce a particle's spin which is parallel to the black-hole 
spin. 

First, for simplicity, we consider a point particle with mass /i that has a spin 
vector = {Si, 0, 0}, and is located at x = {0, 0, R} in the Cartesian coordinates. 
When we discuss perturbations from the spinning particle, we use the following energy- 
momentum tensor. 

rpafi rpafi I rpa/3 . 

mass ~T~ spin ) 

T± = , / , IS = -V, / trSX ^'Z?^ , (31) 
•> v 9 J \ 9 

where we impose a spin supplementary condition (SSC) which determines the center 

of mass of the particle, S^up = 0. Since we focus on the leading order effect of the 

particle's spin, we may consider only the contribution of T^ m . Furthermore, we can 

reduce this energy momentum tensor to the form 

T&n = IS- = -\j dT ^%K) 6W ( X ~ Z ^ ' ^ 

in the leading order of the slow motion approximation. The other components are higher 
order. Here, the spin tensor S l A\ is given by 

sr8 = _^± Si S r<t> = COS COS ^ s ^ = COS0 ^ 

r ' r sin 9 r 

The tensor harmonics coefficients of the stress-energy tensor are calculated as 

A TZ = ^ (sin$9 e y; m (e, $) + ^^ Wje, $)) s(r - r) , 

nsoin i S\ ( ■ ^ r> ^\ cos @ cos $ _ _ _„ _ .\ _ f S(r — R)\ 

lZ = Ju{t+1) ~ ^ m + sin9 °^ (8 ' $) ) ^ ' 
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QfZ = —r4 = - (S%*7J©, $) - cos e cos d^y/je, <&; 



sin 



x* (^) • (34) 

where we set the location of the particle on the z-axis. It should be noted that only 
£ > 1, m — ±1 modes have a non-zero value because of 1^(0,$) ~ (sin0)l m L And 
also, both the even and odd parity modes exist in the metric perturbations from the 
particle's spin. 

We evaluate the I — 2, m — ±1 odd parity perturbations from T^ m in the first 
perturbative order calculation. These perturbations have the leading order effect of the 
particle's spin. The tensor harmonics coefficient of the stress-energy tensor is given by 



(35) 

The wave function is obtained from the Regge- Wheeler equation, 

42 ( i Ad \fr)=S™^(fr),(36) 



d 2 d 2 (r-2M)(r-M) 
dt^ + 'c\* 2 ' Q r 3 



and the source term in the leading order of the slow motion approximation is given by 

5S „» = ^ (~W)Tr S <r - RW) + ^ Hr - m ) • <37) 
where we have extended the interpretation of the particle's location to a sequence of 
quasistatic locations R = R(t). From the above source term, we obtain the Regge- 
Wheeler function 



M« d) (t,r) = ±^S 1 R, (38) 
15 

by using the same Green's function discussed in Subsection 13.21 This wave function has 
the following relation with the black-hole spin effect in Eq. ( 130]) . 

,(l)(odd)/, s Si M (2)(odd) / , s. /qa\ 

i>2±\ (t,r) = -- i/j^l >{t,r). (39) 

3. 6. Evolution of the linear momentum 

Using the results in the previous subsections and the approximated equation of motion in 
the slow motion approximation, d 2 R/dt 2 = —M/R 2 , we can calculate the leading order 
evolution of the linear momentum. The contributions of the gravitational waveforms 
are summarized in Table [TJ The mass, velocity, orbital radius and spin dependence is 
estimated by 

from Eq. ( fit)]) . Also, the direction can be derived from the angular integrations in 
Eq. (HI]). 
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Table 1. Leading order mode contributions to the evolution of the linear momentum 
in the head-on collision of spinning black holes. 



combination 


dependence 


direction 


(£ = 2,m = 0, 


even) 


■ (1 = 


3. 


m 


= 0, even) 


(j?v 7 /R 2 


z 


(£ = 2,m = 0, 


even) 


■ (1 = 


2. 


m 


= ±1, odd) 




x-y 


(£ = 2, m = 0, 


even) 


■(£ = 


2. 


m 


= ±1, odd) 


jj 2 S 2 v 4 /R 4 


x-y 



R - ~ir R • ( 41 ) 



First, we find 

P x = 0. (40) 

This result is also obtained by analyzing the symmetry of the system. In practice, 
we may use the symmetry between m and — m modes in the coefficients of the tensor 
harmonics of the metric perturbation and h^em : and the integration Pf and Pf . 

Next, we discuss the spin independent contribution to the evolution of the linear 
momentum. This contribution arises in the ^-direction as 
■ 16 /x 2 M 2 / - 3 2M 

105 i? 4 \ R 

The above equation is derived from the combination of the first order £ = 2, m = and 
£ = 3, m = even parity modes in Eqs. (fl~8l) and (1201) . 

Finally, for the ^/-direction, we have the leading order contribution of the spin effects 
for the evolution of the linear momentum. 

P y = --li 2 M 2 —J^-?A . (42) 
y 15 P R 5 \M n J K J 

This is calculated from the combination of £ = 2, m = even parity and the £ = 2, 
m = ±1 odd parity modes which include two different wave functions given in Eqs. (130]) 
and 038). 

These results are consistent with the calculation in the post-Newtonian approach [7] . 
Note that Kidder [7] derived the above results for general orbits. On the other hand, 
our calculation discussed here is limited to the head-on collision of spinning black holes 
given in Subsection 13.11 



4. Discussion 



We have discussed gravitational wave recoil effects and analytically derived the leading 
order effects in the evolution of the linear momentum by using the SRWZ formalism. 
This formalism is an extension of the RWZ formalism with a perturbative spin of the 
background black hole. 

In the appendix of [12], we applied the perturbative spin formalism to compute 
the corresponding quasi-normal modes and compare them with those obtained for the 
Kerr black hole for all values of the spin parameter. These results show that the SRWZ 
formalism provide reliable predictions for the spin parameter a/M < 0.3. 
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From the analytic treatment of the gravitational radiation recoil, we confirm the 
leading q 2 (where q = /J>/M) dependence of the large recoils out of the orbital plane |2T| . 
In the black hole perturbation theory, this scaling with mass ratio is trivial, and in [22] 
the spin dependence has also been discussed in detail by solving Teukolsky equation [23] 
numerically. 

The results of the recoil obtained here can be extended to higher PN order in 
the sense of the slow motion approximation where the characteristic orbital velocity 
v <^ 1 (for example, see [21] )• On the other hand, higher order spin effects are 
complicated. In this paper, we have focused only on the couplings between the first order 
perturbations about a Schwarzschild background and the black-hole spin. However, it 
is also necessary to treat the equations of motion with spin. We may use the Teukolsky 
formalism [23] for the spinning large black hole, i.e., use a Kerr background. Taking 
into account the spin of the particle adds a new degree of complication that makes it 
difficult to obtain an analytic expression of the gravitational radiation recoil for general 
orbits and we would need to perform a numerical calculation. 
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